Biomass Valorization through Catalytic Pyrolysis Using Metal-Impregnated Natural Zeolites: From Waste to Resources

Catalytic biomass pyrolysis is one of the most promising routes for obtaining bio-sustainable products that replace petroleum derivatives. This study evaluates the production of aromatic compounds (benzene, toluene, and xylene (BTX)) from the catalytic pyrolysis of lignocellulosic biomass (Pinus radiata (PR) and Eucalyptus globulus (EG)). Chilean natural zeolite (NZ) was used as a catalyst for pyrolysis reactions, which was modified by double ion exchange (H2NZ) and transition metals impregnation (Cu5H2NZ and Ni5H2NZ). The catalysts were characterized by nitrogen adsorption, X-ray diffraction (XRD), ammonium programmed desorption (TPD-NH3), and scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDS). Analytical pyrolysis coupled with gas chromatography/mass spectrometry (Py-GC/MS) allowed us to study the influence of natural and modified zeolite catalysts on BTX production. XRD analysis confirmed the presence of metal oxides (CuO and NiO) in the zeolite framework, and SEM-EDS confirmed successful metal impregnation (6.20% for Cu5H2NZ and 6.97% for Ni5H2NZ). Py-GC/MS revealed a reduction in oxygenated compounds such as esters, ketones, and phenols, along with an increase in aromatic compounds in PR from 2.92% w/w (without catalyst) to 20.89% w/w with Ni5H2NZ at a biomass/catalyst ratio of 1/5, and in EG from 2.69% w/w (without catalyst) to 30.53% w/w with Ni5H2NZ at a biomass/catalyst ratio of 1/2.5. These increases can be attributed to acidic sites within the catalyst pores or on their surface, facilitating deoxygenation reactions such as dehydration, decarboxylation, decarbonylation, aldol condensation, and aromatization. Overall, this study demonstrated that the catalytic biomass pyrolysis process using Chilean natural zeolite modified with double ion exchange and impregnated with transition metals (Cu and Ni) could be highly advantageous for achieving significant conversion of oxygenated compounds into hydrocarbons and, consequently, improving the quality of the condensed pyrolysis vapors.


Introduction
Non-renewable fossil resources such as natural gas, oil, and coal have been the dominant chemical and energy supply for several decades.Global fossil fuel consumption increased from 40,000 terawatt-hours (TWh) in 1965 to 136,000 TWh in 2019 [1].Chemical production also depends on fossil fuels.It is estimated that 12% of the world's oil is Polymers 2024, 16,1912 2 of 18 currently dedicated to producing chemicals that are essential precursors of value-added products such as fertilizers, plastics, rubber, fibers, and solvents [2].Therefore, sustainable alternatives are required to meet the world's energy needs.One of the renewable sources that can fulfill this role is biomass due to its availability, compatibility with current energy infrastructure, and versatility in various applications.
Consequently, numerous studies have explored using agricultural and forestry residues, livestock manure, municipal waste, organic sewage, and industrial refuse as substitutes or oil in their present uses [3].Among these alternatives, lignocellulosic biomass (LCB) is particularly noteworthy, with an annual production rate exceeding 180 billion tons, making it a promising and standout source [2].In 2019, Chile had a forest area of 17 million hectares, of which 14.7 million hectares were native forests and 2.3 million were forest plantations.Pinus radiata and Eucalyptus globulus, with 92.66% of the total, were the most planted species, yielding between 25 and 35 cubic meters per hectare per year [4].Additionally, silvicultural activities and the processing of forest resources generate more than 4 million metric tons of energetically exploitable residues [5].
Biomass can undergo various technical conversion processes to generate different forms of energy.Among these methods, thermochemical conversion techniques such as pyrolysis and hydrothermal liquefaction provide a practical means of producing liquid fuels.These liquid byproducts, called bio-oils, are considered viable alternatives to traditional petroleum fuels for power generation, heating, or extracting valuable chemicals [6].Fast pyrolysis of biomass with a rapid heating rate (>500 • C s −1 ) to intermediate temperatures (400-600 • C) is a promising way to generate bio-oil from the fast decomposition of biomass in the absence of oxygen, by which the short vapor residence time can lead to a high bio-oil yield with few other products like gas and solid char [7].However, the bio-oil produced from pyrolysis is generally an aqueous and highly oxygenated mixture of phenols, acids, and a fraction of hydrocarbons.Furthermore, its high acidity, low stability, corrosive nature, high viscosity, and low vapor pressure are other significant limitations that restrict its direct application in various uses.Therefore, bio-oil upgrading is essential to transform it into a product that can be used and marketed [8].
There are two primary catalytic routes to reduce the oxygen content in bio-oil.The first method, catalytic hydroprocessing, encompasses hydrocracking and hydrodeoxygenation [9].This approach offers the advantage of producing high-quality bio-oil with low oxygen levels and a high H/C ratio.However, it requires a significant hydrogen supply, making the process economically costly [8].The second approach involves catalytic cracking of oxygenated compounds.This process may offer a more cost-effective alternative and could be preferred for upgrading bio-oil [10].The pyrolysis process can be executed in two configurations in catalytic cracking, depending on the catalyst's introduction.The first configuration, termed in situ pyrolysis, entails blending a catalyst with the biomass.Conversely, the second configuration (ex situ pyrolysis) placed the catalyst downstream from the biomass, enabling the resulting pyrolytic vapors to traverse the catalyst bed.Consequently, the catalyst will only come into contact with the volatilized gases [11].In in situ pyrolysis, the generated pyrolytic vapors cannot come into contact with a significant amount of catalysts, thus requiring a higher biomass/catalyst ratio to achieve better deoxygenation activity [12].In the catalytic pyrolysis of hybrid poplar, Wang et al. [13] used an HZSM-5 zeolite as a catalyst, finding that aromatics formed from ex situ pyrolysis were predominantly monocyclic aromatics, whereas in situ pyrolysis favored polycyclic aromatics.
Using zeolite-based catalysts decomposes and dehydrates biomass polymers due to their high selectivity towards producing aromatic compounds and their ability to limit the formation of unwanted products, such as carboxylic acids and oxygenated compounds [14].In this context, renewable aromatic compounds are crucial foundational components capable of partially substituting petrochemical resources in producing biofuels and chemicals.They have the potential to diminish reliance on oil and mitigate environmental repercussions by emitting fewer SOx and CO 2 emissions, thus fostering the growth of sustainable economies [15].Synthetic zeolites (Y, ZSM-5, and H-ZSM-5) are the preferred catalysts in pyrolysis oil upgrading processes, but they have a limitation due to their high cost (e.g., ZSM-5 → 120 USD•kg −1 , Y → 180 USD•kg −1 ).In contrast, natural zeolites are inexpensive and abundant; for example, natural clinoptilolite has an average value of 242 USD• metric ton −1 [16].Besides the lower cost of extraction, natural zeolites have some advantages over synthetic zeolites, like reduced environmental impact in their extraction, as synthetic zeolite production involves chemicals and intensive industrial processes.Additionally, the production of natural zeolites generally requires less energy than zeolite synthesis, which can result in a lower energy footprint.Finally, natural zeolites require less post-processing for industrial applications, simplifying preparation procedures and reducing associated costs [17,18].
The drawback of natural zeolites is their lower performance in pyrolysis compared to synthetic ones, which can be significantly improved through appropriate treatments [19].The ionic exchange procedure enhances Brönsted acid sites (bridging hydroxyls connected to aluminum atoms in the framework) within the internal structure of zeolites, leading to bio-oil deoxygenation [20].Additionally, the impregnation of metals and metal oxides into zeolites reduces the formation of acids and oxygenated compounds in catalytic pyrolysis since the metal particles can provide sites for the formation and breaking of C-C bonds and the acid sites of the zeolites can act in the isomerization of olefins [21].In a study conducted by Veses et al. [22], the performance of various metal cations/ZSM-5 catalysts in bio-oil upgrading through ex situ pyrolysis of biomass was compared.They demonstrated that Ni/ZSM-5 and Cu/ZSM-5 catalysts outperformed Mg/ZSM-5 and H-ZSM-5 catalysts.Expressly, the results indicated that Ni/ZSM-5 and Cu/ZSM-5 yielded bio-oil containing approximately 35% and 31% hydrocarbons, respectively, while Mg/ZSM-5 and H-ZSM-5 catalysts produced bio-oil samples with 29% hydrocarbons.Rajić et al. [23] compared the catalytic activity of clinoptilolite containing NiO, ZnO, and Cu 2 O nanoparticles in the pyrolysis of hardwood lignin.They found that the production of phenolic compounds primarily depended on the metal cations in the zeolite and not on the type of acid sites present.The bio-oil produced contained phenols in the range of 39% for clinoptilolite with ZnO, 50% for clinoptilolite with Cu 2 O, and 54% for clinoptilolite with NiO.They attribute the higher yield of phenols in hardwood pyrolysis to using clinoptilolite with NiO as a catalyst for the interaction between the metal nanoparticles and the zeolite during the reaction.Regardless of the prominent results obtained, the significant role of the catalyst in pyrolysis product distribution is inconclusive and requires further analysis to elucidate the reaction pathways induced by the catalyst addition (clinoptilolite).
Despite the extensive literature on the effect of synthetic zeolite-based catalysts in lignocellulosic biomass pyrolysis, there is currently no consensus on the influence of natural zeolite application as their replacement.Previous evidence indicates that the formation of BTX is promoted due to the synergy between lignocellulosic biomass components and zeolite-based catalysts.Therefore, this study aims to evaluate the composition of compounds formed during the catalytic pyrolysis of Pinus radiata and Eucalyptus globulus from the Bio-Bio Region, Chile, using Chilean natural zeolite impregnated with transition metals (Cu and Ni) and their influence on the formation of aromatic compounds: benzene, toluene, and xylene (BTX), as well as the reduction in oxygenated compounds.

Biomass
Investigaciones Forestales Bioforest S.A., based in Bío Bío, Chile, supplied biomass samples (PR and EG) as wood chips.These chips were composed entirely of wood without bark, and the specific location of the sample within the tree trunk (sapwood/heartwood) was not considered.The samples were ground and sieved to attain an average particle size of 0.33 mm.Following a previously documented procedure, they were stored in a desiccator until further use [24].Chilean natural zeolite (53% clinoptilolite, 40% mordenite, and 7% quartz) was provided by "Minera Formas".The zeolite sample was ground and sieved with an average particle size of 0.3 mm, and then rinsed with ultrapure water, oven-dried at 125 • C for 24 h, and finally stored in a desiccator until further use [25].

Catalyst Modification
The Chilean natural zeolite (NZ) was modified by ion exchange using ammonium sulphate [(NH 4 ) 2 SO 4 ] (0.1 mol•dm −3 ), according to the previous work reported by Alejandro et al. [26].The exchanged zeolite was denoted as (H2NZ).Finally, the sample was dried in an oven at 125 • C for 24 h and subjected to thermal degassing with a flow of N 2 inert gas (100 cm 3 •min −1 ) for two hours with a heating rate of 10 • C•min −1 .The H2NZ sample was impregnated with transition metals (Cu and Ni) using metal salts precursors: nickel nitrate Ni(NO 3 ) 2 •6H 2 O and copper nitrate Cu(NO3) 2 •3H 2 O, through the incipient wet impregnation technique, as reported previously [27] and considering a metallic salts concentration of 0.05 mol•dm −3 [28,29].Finally, the samples were calcined at 500 • C for 4 h and then cooled until reaching room temperature.Such samples were named Cu5H2NZ and Ni5H2NZ, respectively.

Biomass Characterization
The elemental analysis of PR and EG samples followed the ASTM D5373 Standard method [30] and was conducted using a Leco CHNS 628 elemental analyzer (St.Joseph, MI, USA).Proximate analyses, which included determining moisture content, volatile matter, ash content, and fixed carbon, were carried out in a muffle furnace, adhering to the ASTM D3172 Standard method [31].The oxygen (O) content was calculated by the difference between carbon (C), hydrogen (H), and nitrogen (N) measurements.
Around 300-400 mg of PR and EG with 45-60 mesh were considered to quantify extractives, using a Soxhlet system with an acetone/water solution (9/1) for 16 h, following the TAPPI T280 pm-99 method [32], under the procedure described by Aguayo et al. [33].The methods for determining holocellulose, α-cellulose, and lignin can be found in previous research [34].

Physico-Chemical Characterization of Natural and Modified Zeolite Samples
Both natural and modified zeolites were characterized using various physico-chemical techniques.X-ray powder diffraction (XRD) analysis was used to confirm the modifications of the zeolite structure.X-ray diffraction patterns were obtained using a Bruker D8 Advance equipment (Billerica, MA, USA) operated at 20 mA and 40 kV with a copper cathode lamp (λ = 1.541Å).XRD patterns were collected in the angular range 6-70 • , a step size of 0.02 and a time/step of 0.2 s.The Scherrer equation (Equation ( 1)) [35] can be estimated from the width of the X-ray diffraction peaks where D = Average crystal size, λ = X-ray wavelength (0.154 nm), θ = Bragg angle, and B = FWHM in radians.
In order to compare the crystallinity from the natural and modified zeolites, the relative crystallinity was determined according to ASTM D5357 [36] recommendations (Equation ( 2)).This value compares the peak intensity corresponding to the planes of the crystalline zeolite with the modified zeolite peaks [37].
where crystallinity = relative crystallinity (%), S sample = sum of integral peak intensities for the sample, and S reference = sum of integral peak intensities for the reference.The surface morphology of the natural and modified zeolites was examined by SEM-JEOL JSM-IT300 scanning electron microscopy (JEOL, Tokyo, Japan) at 5-20 kV and 50 Pa.Energy-dispersed X-ray spectroscopy (EDS) was performed to confirm the presence of the impregnated transition metals on the surface.Specific surface areas were determined by nitrogen adsorption and desorption at 77 K using a Micromeritics Gemini 2370 apparatus (Norcross, GA, USA), using 0.57 g of samples previously outgassed under nitrogen flow at 623 K for 2 h.The total area was calculated using the Brunauer-Emmett-Teller (BET) model with the Quantachrome, NOVATOUCH LX1 model.The total pore volume was obtained at P/P 0 = 0.99.Micropore area and micropore volume were calculated by the t-plot method using the adsorption data range of 0.2 < P/P 0 < 0.6.The pore size distribution was derived from the adsorption branch data of each isotherm using the Barrett-Joyner-Halenda (BJH) method [38].
The TGA experiments were carried out in a TGA Q500 thermogravimetric analyzer (TA Instruments, New Castle, DE, USA) to evaluate the thermal stability of the natural and modified zeolite.Approximately 10 ± 1 mg of the sample was placed on the thermobalance and heated from room temperature to 650 • C at a heating rate of 10 • C•min −1 .The experiments were conducted at atmospheric pressure with an N2 flow (50 mL•min −1 ).
Ammonia thermal programmed desorption (NH 3 -TPD) contributes to the investigation of acid sites of the zeolite samples [39].Samples (0.055 g) were degassed under nitrogen flow (0.05 L•min −1 ), heated to 550 • C at 10 • C•min −1 , and finally cooled down to 125 • C to saturated with ammonia.After the samples were heated up to 550 ), the TPD system used a thermal conductivity detector (TCD) to determine the desorbed ammonia.The changes in conductivity vs. temperature were recorded since the temperature ranges in which NH 3 desorbs are directly related to the strength of the acid sites of the zeolite, and the amount of chemisorbed NH 3 is proportional to the number of acid sites per unit mass of the adsorbent [40].

Catalytic Fast Pyrolysis of Biomass
The catalytic biomass pyrolysis was carried out in an analytical micro-pyrolysis reactor (Pyroprobe 5200HPR, CDS Analytical Co., Ltd., Oxford, PA, USA) connected in line with a gas chromatograph equipped with a mass spectrometer system, GC/MS (Clarus 690, QS8.Perkin Elmer, Waltham, MA, USA).The pyrolysis reactor consisted of a quartz tube heated by a platinum filament.Approximately 0.5 ± 0.1 mg of each sample was weighed by a microbalance (AD 6000 Ultra MicroBalance Perkin Elmer, Waltham, MA, USA).According to previous analyses, the experiments were conducted at 550 • C to guarantee maximum biomass decomposition and to achieve the highest pyrolysis oil yield [41].He (pure 99.996%, BOC) was used as a carrier gas.The compounds in the pyrolysis-evolved gas were identified by referencing the NIST 2017 library and TurboMass 6.1.0software.The ex situ configuration was chosen in the catalytic pyrolysis experiments [42].In this configuration, the catalyst has no contact with the biomass, which allows for a more significant interaction between the pyrolysis vapors and the active sites of the catalysts.Additionally, Table 1 shows the conditions of catalytic biomass pyrolysis.A semi-quantitative approach was employed to estimate the concentration of aromatic compounds, utilizing the ratio of absolute peak area to sample mass (as depicted in Equation ( 3)) [43].
The selective yield of individual BTX compounds (benzene, toluene, and xylene) within the total aromatic compounds was determined using Equation (4) [44].

Biomass Characterization
Table 2 presents the raw materials' proximate analysis, ultimate analysis, and chemical composition of biomass samples (PR and EG).The proximate analysis revealed no significant differences among these wood species, indicating a high volatile content that supports their suitability for pyrolysis.Notably, the ash content in these samples was below 1%, eliminating concerns about the potential catalytic effects of ashes during pyrolysis [45].The ultimate analysis poses the PR and EG at the same level as other woody biomass previously reported in the literature [46].The high carbon (48.02-47.76%)and low nitrogen (0.29-0.09%) contents further affirm the appropriateness of PR and EG for thermochemical conversion processes and suggest that this raw material is unlikely to generate NOx and SOx emissions.

Crystallinity
Figure 1 shows the X-ray diffraction patterns for metal (Cu, Ni) impregnated zeolites.It can be revealed in the XRD that the natural zeolite sample is highly crystalline, showing characteristic peaks of mordenite (M), clinoptilolite (C), and quartz (Q), which have been previously reported [23].The characteristic peaks of quartz are observed at (21.0 • and 26.7 • ), just as reported by Gurevich et al. [20].It can be seen that the modifications made (double ion exchange and impregnation of transition metals) did not affect the crystallinity of the catalysts, as already reported by Trisunaryanti et al. [48].
International Centre for Diffraction Data (ICDD) reference codes 00-045-0937 and 01-7390.The metal oxides in the metal-modified zeolites were mainly dispersed on the face of H2NZ and partially anchored in the pores of the catalysts, thus preserving structure [50].Thus, the low content of Cu and Ni doping reduced this agglomeration the zeolite surface [51].Using the Scherrer equation, crystal size can be estimated from the width of X diffraction peaks.The results of crystal size measurements for each impregnated m oxide phase are described in Table 3.The increase in crystal size is reflected in the decre in pore size and surface area of the catalyst due to transition metal impregnation [35], this is possible since the more impregnated metal into the pores of the carrier, the m the pores of the carrier with smaller fingers will be blocked [35].
Table 3 also shows the relative crystallinity of zeolites.The relative crystallinity of natural and modified zeolite samples was calculated, confirming that the highest inten was obtained in H2NZ, which was taken as the reference.During ion exchange, the cat present in the channels and cavities of the zeolite are replaced.This process can rearra the crystalline structure, eliminating defects and improving the organization of the cry lattice.This rearrangement can result in a more ordered structure and, consequently increase in the crystallinity.This might explain why the crystallinity of the ion-exchan modified zeolite is higher than that of the natural zeolite [52].The ion exchange pro can also help eliminate impurities in the natural zeolite.These impurities may have cau distortions or defects in the original crystalline structure.By removing these impuri the structure can become more homogeneous and ordered, which is reflected in a hig The appearance of new peaks in the Cu5H2NZ and Ni5H2NZ diffractograms indicates the presence of the oxide metals on the zeolite surface and confirms the results obtained by SEM-EDS. Figure 1 shows a CuO phase at peaks of 35.10 • , 39.30, 53.44 • , and 61.40, which can be indexed to (002), ( 200), (020), and (−311) planes of crystalline CuO [8].The peaks identified for nickel oxide (NiO) were found at positions 2θ: 37.25, 43.29, and 62.85 [49], which were indexed to (111), (200), and (220) planes of crystalline NiO.These results were consistent with the standard values of CuO and NiO and well-matched with International Centre for Diffraction Data (ICDD) reference codes 00-045-0937 and 01-089-7390.The metal oxides in the metal-modified zeolites were mainly dispersed on the surface of H2NZ and partially anchored in the pores of the catalysts, thus preserving the structure [50].Thus, the low content of Cu and Ni doping reduced this agglomeration on the zeolite surface [51].
Using the Scherrer equation, crystal size can be estimated from the width of X-ray diffraction peaks.The results of crystal size measurements for each impregnated metal oxide phase are described in Table 3.The increase in crystal size is reflected in the decrease in pore size and surface area of the catalyst due to transition metal impregnation [35], and this is possible since the more impregnated metal into the pores of the carrier, the more the pores of the carrier with smaller fingers will be blocked [35].Table 3 also shows the relative crystallinity of zeolites.The relative crystallinity of the natural and modified zeolite samples was calculated, confirming that the highest intensity was obtained in H2NZ, which was taken as the reference.During ion exchange, the cations present in the channels and cavities of the zeolite are replaced.This process can rearrange the crystalline structure, eliminating defects and improving the organization of the crystal lattice.This rearrangement can result in a more ordered structure and, consequently, an increase in the crystallinity.This might explain why the crystallinity of the ion-exchanged modified zeolite is higher than that of the natural zeolite [52].The ion exchange process can also help eliminate impurities in the natural zeolite.These impurities may have caused distortions or defects in the original crystalline structure.By removing these impurities, the structure can become more homogeneous and ordered, which is reflected in a higher crystallinity index.Finally, ion exchange can improve the mobility of ions within the zeolite structure.Better ionic mobility can facilitate a more efficient rearrangement of atoms in the crystal lattice, promoting higher crystallinity [53][54][55].The results demonstrate that the transition metals impregnation in zeolites (Cu5H2NZ and Ni5H2NZ) slightly decreased crystallinity compared to parent zeolites.However, this reduction did not damage the crystalline structure of the zeolite [50].

Chemical Composition and Surface Morphology
Energy-dispersive X-ray spectroscopy (EDS) was conducted in order to identify the compensating cations (Na, K, Mg, Ca, Mg, Fe, and Ti) in the zeolite framework, as well as the impregnation of metals (Cu and Ni), and the results are shown in Table 4.The catalyst underwent an ion exchange process with NH 4 + , decreasing Na + , K + , and Ca 2+ concentrations.Similar results were obtained in a previous work reported by Alejandro et al. [26].These cations were found to have a lower affinity for the zeolite framework than NH 4 + .However, the Mg 2+ cation was hardly affected by the ion exchange and had a higher affinity for clinoptilolite than NH 4 + [56].Generally, the surface concentration of compensation cations decreased due to replacing most of the alkaline and alkaline earth cations with NH 4 + , which later decomposed to H + to generate Brönsted acid sites.Na ions in natural clinoptilolite zeolites are weakly bonded to the zeolite framework and are easy to remove, as confirmed by Ates and Hardacre [55].During chemical modification, calcium cations move into the solution.In contrast, magnesium and potassium cations are less exchanged, probably due to the difficulty in accessing magnesium and potassium cations because of their larger radius compared to Ca 2+ cations (for instance, Mg 2+ has an atomic radius of 0.16 nm, while K + and Ca 2+ have atomic radii of 0.133 nm and 0.104 nm, respectively) [57].The results of the copper and nickel-loaded samples show that metal impregnation indeed occurred and replaced only the ions of the compensation cations in the clinoptilolite phase.Similar results were observed by Rajić et al. [23].Figure 2 shows the scanning electron microscope (SEM) images before and after impregnation with 2000X magnification.The natural zeolite (NZ) (Figure 2a) and natural with double ion exchange (H2NZ) (Figure 2b) present an irregular surface morphology without presenting significant changes between the zeolites.Figure 2c,d shows SEM images of the zeolites impregnated with transition metals Cu and Ni, respectively.The presence of white dots on the surface is noted, which indicates the presence of metals in the zeolite, and the surfaces become more regular compared to NZ and H2NZ, which may be due to the level of dispersion and uniform distribution of the metals on the surface of the zeolite.Similar observations were reported by Trisunaryanti et al. [58].

Textural Properties
Textural properties of natural and modified zeolite samples are listed in Table 5. Surface areas (SBET) were calculated by applying the BET adsorption model to nitrogen adsorption data.The results indicate that an ion exchange process with NH4 + decreased the surface area without significant Si/Al ratio changes.These results are consistent with those published by Alejandro et al. [59].This decrease in the surface area of Cu5H2NZ and Ni5H2NZ can be attributed to the incorporation of metal ions onto the surface of zeolites or inside the pores of the zeolite structure.This incorporation of transition metals could potentially block the zeolite pores, as indicated by Veses et al. [60].
Furthermore, the average pore size decreased for this catalyst, implying that more micropores were created in the Cu and Ni impregnation [20].Metal impregnation (Cu, Ni) onto the zeolite framework decreased its surface area.These phenomena may be caused by metal distribution on the catalyst surface.Some of the metal atoms may block the pore mouth of the zeolite [48].Furthermore, it was observed that the total pore volume in Cu5H2NZ and Ni5H2NZ was lower compared to the NZ.This decrease in pore volume due to metal loading suggests pore blockage, which can be attributed to the effective deposition of metals on the zeolite surface or their potential presence within the zeolite pores.

Textural Properties
Textural properties of natural and modified zeolite samples are listed in Table 5. Surface areas (SBET) were calculated by applying the BET adsorption model to nitrogen adsorption data.The results indicate that an ion exchange process with NH 4 + decreased the surface area without significant Si/Al ratio changes.These results are consistent with those published by Alejandro et al. [59].This decrease in the surface area of Cu5H2NZ and Ni5H2NZ can be attributed to the incorporation of metal ions onto the surface of zeolites or inside the pores of the zeolite structure.This incorporation of transition metals could potentially block the zeolite pores, as indicated by Veses et al. [60].

Sample
S BET (m 2 •g −1 ) Pore Size (nm) V micropore (cm 3 •g −1 ) V mesopore (cm 3  Furthermore, the average pore size decreased for this catalyst, implying that more micropores were created in the Cu and Ni impregnation [20].Metal impregnation (Cu, Ni) onto the zeolite framework decreased its surface area.These phenomena may be caused by metal distribution on the catalyst surface.Some of the metal atoms may block the pore mouth of the zeolite [48].Furthermore, it was observed that the total pore volume in Cu5H2NZ and Ni5H2NZ was lower compared to the NZ.This decrease in pore volume due to metal loading suggests pore blockage, which can be attributed to the effective deposition of metals on the zeolite surface or their potential presence within the zeolite pores.

Thermal Analysis
Natural zeolites impregnated with Cu and Ni must be thermally and chemically stable to be considered as potential catalysts.Figure 3      The NH 3 desorption occurs across a wide temperature range, indicating weak and strong acid sites in all the samples.The peaks with a maximum below 325 • C correspond to weak acidic sites, as those above 325 • C represent strong acidic sites [63].An increase in the total acidic sites was observed for impregnated zeolites.It can be associated with metals (Cu and Ni) that are highly coordinated (they can bond with the ligand due to the arrangement of the electrons in the unoccupied d orbital), which has low-energy orbitals and can act as a weak acidic site [64].
Table 6 shows the strength of the acid sites in the natural and modified zeolite samples.A high number and density of acid sites are desirable for deoxygenation reactions and the production of aromatic compounds during the pyrolysis of lignocellulosic biomass [65].H2NZ posed stronger acid sites than NZ, promoted by Brönsted acid sites associated with hydroxyl groups, generated after ion exchange with NH 4 + and thermal treatment [12].Impregnation introduced Cu 2+ and Ni 2+ cations, present as CuO and NiO oxides, as confirmed in the XRD analysis.The presence of these oxides led to an increase in strong acid sites at lower temperatures, which aligns with what was published by Veses et al. [22].A GC/MS analysis was employed to identify obtained compounds during the catalytic pyrolysis of PR and EG using NZ, H2NZ, Cu5H2NZ, and Ni5H2NZ as catalysts for the biomass-catalyst combinations: 1/1, 1/2.5, and 1/5.The family compounds detected are shown in Table S1 of the Supplementary Information.The compounds and compound families were quantified using a semi-quantitative approach that determined the relative area percentage of all identified constituents, which were categorized into different groups (phenols, acids, aldehydes, furans, ketones, aromatics, and esters) based on previous reports [66].Non-catalytic pyrolysis assays were also conducted to compare catalyst performance.
Oxygenated compounds, such as acids, aldehydes, esters, and ketones, are considered undesirable fractions for energy production, while hydrocarbons and alcohols are regarded as desirable products for biofuel production.Phenols and furans are also seen as highvalue chemicals [67].The data indicate that non-catalytic pyrolysis led to the production of condensate vapors with a high content of oxygenated compounds, primarily yielding approximately 5.40% aldehydes, 8.74% acids, 12.09% esters, and 18.75% ketones in PR and 14.51% aldehydes, 15.31% acids, 6.37% esters, and 13.78% ketones in EG.Other oxygenated compounds were also present, including 5.72% alcohols, 4.62% furans, and 32.77% phenols in PR.On the other hand, 4.61% alcohols, 4.94% furans, and 24.06% phenols were observed when EG was used.This higher proportion of phenols was expected to result primarily from the pyrolysis of the lignin component of the biomass.At the same time, the production of acids and ketones could be chiefly attributed to the thermal degradation of cellulose, hemicellulose, and some portion of lignin [68].
Further analysis using catalyst samples (NZ, H2NZ, Cu5H2NZ, and Ni5H2NZ) revealed a significant reduction in the formation of acids, esters, ketones, alcohols, and phenols [69].There was a pronounced reduction of ketones for PR and less evident for EG.Likewise, Cu5H2NZ samples promoted 22% phenols and 10% ketones, and Ni5H2NZ 27.5% phenols and 12.5% ketones in EG pyrolysis.Such results suggest that the catalysts favored deoxygenation pathways such as dehydration, decarboxylation, and decarbonylation to convert oxygenated compounds [70].The presence of acidic compounds can be attributed to cellulose and hemicellulose depolymerization due to dehydration, C-C bond cleavage, and direct C-O bond cleavage occurring at the acidic sites of catalysts [68].
Thus, the increased deoxygenation activity of Ni5H2NZ, compared to Cu5H2NZ, can be attributed to its higher surface area, porosity, and acid site density.A predominant deoxygenation reaction for the Cu5H2NZ sample could be possibly associated with acid decarbonylation to form aldehydes at copper cations.On the other hand, the preferred mechanism for the O-removal using Ni5H2NZ samples seems to involve decarbonylation and decarboxylation reactions occurring at the acid sites formed after metal integration.Although a dominant reaction mechanism could not be assured, the lower incorporation of cations at ion exchange positions could potentially assess the reduced deoxygenation rate [60].Additionally, they demonstrated that the increased aromatics fraction correlates with the increased number of accessible acid sites at the mesopore surface, where preferential occurrence of decarbonylation reactions is essential for enhancing aromatic production.The deoxygenation of condensed vapors increases with the number of weak acid sites created by incorporating cations on the external surface and mesoporous walls, confirming previously reported results by Veses et al. [60], who demonstrated a direct relationship between the deoxygenation rate and the weak acid sites increment [11].
The higher production of aromatic compounds was reported for EG (around 30%), while in PR, aromatics were close to 20%.Additionally, aromatics production is higher at a B/C ratio of 1/5 in PR and 1/2.5 in EG for all catalysts.The deoxygenation and increased proportion of aromatics are related to a higher catalyst amount, allowing more active sites to interact with pyrolytic vapors, resulting in greater deoxygenation.
Figure 5 displays the compound families formed during biomass pyrolysis and catalytic biomass pyrolysis.The single color bars correspond to the 1/1 biomass/catalyst, the bars with two colors (colored and white) correspond to 1/2.5, and the bars' horizontal line corresponds to the 1/5 condition.The black bars correspond to the trials without a catalyst (WOC), the red bars correspond to NZ, the blue bars to H2NZ, the green bars to Cu5H2NZ, and the brown bars to Ni5H2NZ.The components detected in the condensed vapors were divided into three groups: oxygenates, phenols, and aromatics.Aromatics were divided into monocyclic aromatic hydrocarbons (MAHs), mainly benzene, toluene, and xylene (BTX), and polycyclic aromatic hydrocarbons (PAHs).Moreover, oxygenate production was based on alcohols, aldehydes, carboxylic acids, esters, furans, and ketones.27.5% phenols and 12.5% ketones in EG pyrolysis.Such results suggest that the catalysts favored deoxygenation pathways such as dehydration, decarboxylation, and decarbonylation to convert oxygenated compounds [70].The presence of acidic compounds can be attributed to cellulose and hemicellulose depolymerization due to dehydration, C-C bond cleavage, and direct C-O bond cleavage occurring at the acidic sites of catalysts [68].Thus, the increased deoxygenation activity of Ni5H2NZ, compared to Cu5H2NZ, can be attributed to its higher surface area, porosity, and acid site density.A predominant deoxygenation reaction for the Cu5H2NZ sample could be possibly associated with acid decarbonylation to form aldehydes at copper cations.On the other hand, the preferred mechanism for the O-removal using Ni5H2NZ samples seems to involve decarbonylation and decarboxylation reactions occurring at the acid sites formed after metal integration.Although a dominant reaction mechanism could not be assured, the lower incorporation of cations at ion exchange positions could potentially assess the reduced deoxygenation rate [60].Additionally, they demonstrated that the increased aromatics fraction correlates with the increased number of accessible acid sites at the mesopore surface, where preferential occurrence of decarbonylation reactions is essential for enhancing aromatic production.The deoxygenation of condensed vapors increases with the number of weak acid sites created by incorporating cations on the external surface and mesoporous walls, confirming previously reported results by Veses et al. [60], who demonstrated a direct relationship between the deoxygenation rate and the weak acid sites increment [11].
The higher production of aromatic compounds was reported for EG (around 30%), while in PR, aromatics were close to 20%.Additionally, aromatics production is higher at a B/C ratio of 1/5 in PR and 1/2.5 in EG for all catalysts.The deoxygenation and increased proportion of aromatics are related to a higher catalyst amount, allowing more active sites to interact with pyrolytic vapors, resulting in greater deoxygenation.
Figure 5 displays the compound families formed during biomass pyrolysis and catalytic biomass pyrolysis.The single color bars correspond to the 1/1 biomass/catalyst, the bars with two colors (colored and white) correspond to 1/2.5, and the bars' horizontal line corresponds to the 1/5 condition.The black bars correspond to the trials without a catalyst (WOC), the red bars correspond to NZ, the blue bars to H2NZ, the green bars to Cu5H2NZ, and the brown bars to Ni5H2NZ.The components detected in the condensed vapors were divided into three groups: oxygenates, phenols, and aromatics.Aromatics were divided into monocyclic aromatic hydrocarbons (MAHs), mainly benzene, toluene, and xylene (BTX), and polycyclic aromatic hydrocarbons (PAHs).Moreover, oxygenate production was based on alcohols, aldehydes, carboxylic acids, esters, furans, and ketones.
The analytical results indicate that the presence of catalysts (NZ, H2NZ, Cu5 H2NZ, and Ni5H2NZ) can enhance aromatics production from biomass catalytic conversion.Ni5H2NZ can better promote the formation of BTX and reduce the formation of oxygenated compounds than Cu5H2NZ.
During biomass catalytic pyrolysis, oxygenated pyrolytic vapors with large molecules were first broken down into smaller hydrocarbon precursors (including monofunctional furans, ketones, and light phenols) by the metallic oxides, which favored both diffusion and catalytic reactions in the zeolite, resulting in the improved formation of aromatic compounds [74].
The decomposition of cellulose and hemicellulose generates oxygenated compounds, including anhydrosugars, acids, and furfural.These compounds undergo decarboxylation, ketonization, and aldol condensation reactions over the metallic oxides to produce furans and ketones, which are then transformed into aromatic compounds [75].Furans are indispensable intermediates for the generation of BTX, which the catalyst promotes.
The "hydrocarbon pool", a mass of olefins, is produced through oligomerization and can be converted into aromatics via aromatization over the catalyst.The generated olefins can also react with different furans to form aromatics through the Diels-Alder reaction [8].
Ketones also undergo aromatization in the catalyst channels through self-condensation, decarboxylation, and decarbonylation reactions.Moreover, light oxygenates could be deoxygenated and cracked to C2 to C6 olefins over the catalyst, which would be aromatized into benzene and other aromatics [76].

Aromatics Compounds
The highest production of aromatic compounds was observed for Pinus radiata (20.89% with Ni5H2NZ, biomass/catalysts = 1/5) and Eucalyptus globulus (35.53% with Cu5H2NZ, biomass/catalysts = 1/2.5),with a generally higher production of aromatics in the catalytic pyrolysis of EG than PR.Then, aromatics production was at its maximum at a B/C ratio of 1/2.5 for PR and 1/5 for EG.Deoxygenation and aromatics increments are related to a higher density of active sites in the catalysts interacting with pyrolytic vapors.Additionally, the increase in the aromatic fraction could be related to more available strong acid sites on the catalysts' surface, which are responsible for decarbonylation reactions and essential for improving aromatic production.Deoxygenation increased with the density of weak acid sites formed by incorporating cations on the external surface and mesoporous walls, as reported elsewhere [11,60].
During catalytic pyrolysis, a biomass/catalyst ratio of 2.5 proved the most significant in achieving maximum aromatic formation.A higher amount of catalyst contributes to more active sites involved with pyrolytic vapors, resulting in deoxygenation of condensable vapors and enhanced aromatic production compared to biomass/catalyst ratios of 1/1 and 1/5.All catalysts showed a reduction in the proportion of oxygenated compounds.Adding Cu or Ni also improved the efficiency of the catalysts in deoxygenating condensable vapors.The proportion of acidic compounds was considerably reduced.These results align with those published by Tian et al. [44].During catalytic pyrolysis, a biomass/catalyst ratio of 2.5 proved the most significant in achieving maximum aromatic formation.A higher amount of catalyst contributes to more active sites involved with pyrolytic vapors, resulting in deoxygenation of condensable vapors and enhanced aromatic production compared to biomass/catalyst ratios of 1/1 and 1/5.All catalysts showed a reduction in the proportion of oxygenated compounds.Adding Cu or Ni also improved the efficiency of the catalysts in deoxygenating condensable vapors.The proportion of acidic compounds was considerably reduced.These results align with those published by Tian et al. [44].
This efficient aromatic production can be attributed to the excellent catalytic activity of Cu and Ni cations and physicochemical properties such as higher surface area and an increased number of acidic sites.Additionally, this production can be attributed to the aromatization activity carried out by the acidic sites present in the catalysts and other deoxygenation pathways, such as dehydration, decarboxylation, and decarbonylation, which were efficiently catalyzed by the catalysts [8].
Table S2 shows the effect of the catalysts on BTX yields, and Table S3 presents the effect of the catalysts on the aromatic's selectivity.The total BTX selectivity increased after adding the catalysts.These yields obtained with (NZ, H2NZ, Cu5H2NZ, and Ni5H2NZ) were attributed to their suitable pore structure and acidic strength.Aromatics were derived from the selective catalytic fast pyrolysis of oxygenates (especially G and S-type phenols) through dihydroxylation and demethoxylation reactions [80].These intermediate oxygenates entered the zeolite channels and underwent a series of deoxygenation reactions at the acid sites [43].Among the BTX compounds, toluene was the most prevalent for Pinus radiata and Eucalyptus globulus.Thus, the selective yields of toluene and xylene in the 1/2.5 ratio were higher than those of benzene because toluene and xylene can be directly produced from the dehydration (-OH) and demethoxylation (-OCH3) of oxygenated precursors (especially phenolics).In contrast, an additional step is required to produce benzene compared to the steps needed for toluene and xylene: the demethylation (-CH3) of toluene and xylene [81].Therefore, obtaining benzene is more challenging than toluene and xylene, hence the lower yield [80].This efficient aromatic production can be attributed to the excellent catalytic activity of Cu and Ni cations and physicochemical properties such as higher surface area and an increased number of acidic sites.Additionally, this production can be attributed to the aromatization activity carried out by the acidic sites present in the catalysts and other deoxygenation pathways, such as dehydration, decarboxylation, and decarbonylation, which were efficiently catalyzed by the catalysts [8].
Table S2 shows the effect of the catalysts on BTX yields, and Table S3 presents the effect of the catalysts on the aromatic's selectivity.The total BTX selectivity increased after adding the catalysts.These yields obtained with (NZ, H2NZ, Cu5H2NZ, and Ni5H2NZ) were attributed to their suitable pore structure and acidic strength.Aromatics were derived from the selective catalytic fast pyrolysis of oxygenates (especially G and S-type phenols) through dihydroxylation and demethoxylation reactions [80].These intermediate oxygenates entered the zeolite channels and underwent a series of deoxygenation reactions at the acid sites [43].Among the BTX compounds, toluene was the most prevalent for Pinus radiata and Eucalyptus globulus.Thus, the selective yields of toluene and xylene in the 1/2.5 ratio were higher than those of benzene because toluene and xylene can be directly produced from the dehydration (-OH) and demethoxylation (-OCH 3 ) of oxygenated precursors (especially phenolics).In contrast, an additional step is required to produce benzene compared to the steps needed for toluene and xylene: the demethylation (-CH 3 ) of toluene and xylene [81].Therefore, obtaining benzene is more challenging than toluene and xylene, hence the lower yield [80].

Conclusions
This study evaluated the influence of metal-supported Chilean natural zeolites in the analytical pyrolysis of Pinus radiata and Eucalyptus globulus.Natural zeolite was modified to incorporate transition metals (Cu and Ni) on the surface through incipient wet impregnation.SEM-EDS analyses confirmed the successful impregnation (6.20% for Cu5H2NZ and 6.97% for Ni5H2NZ).XRD analysis showed the zeolite framework's incorporated metal oxides (CuO and NiO).

Figure 3 .
Figure 3. Thermal analysis of zeolite samples (a) TG and (b) DTG.3.2.5.Acid Site StrengthNH3-TPD characterizes the strength of accessible acid sites in the catalysts, as Figure4represents.

Table 2 .
Characterization results of biomass samples.

Table 3 .
Crystal size in zeolite samples.

Table 4 .
Physico-chemical characterization of the Chilean natural zeolite and impregnated.
nd: Not detected.

Table 6 .
Physico-chemical characterization of the natural and impregnated zeolites.